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Abstract. The equation of state of cold baryonic matter is studied within a relativistic 
mean-field model with hadron masses and coupling constants depending on a scalar field. We 
demonstrate that if the effective nucleon mass stops to decrease with a density increase at 
densities n > n* > no, where no is the nuclear saturation density, the equation of state stiffens 
for these densities and the limiting neutron star mass increases. The stabilization of the nucleon 
mass can be realised if in the equation of motion for the scalar mean-field there appear a term 
sharply varying in a narrow vicinity of the field value corresponding to the density n*. We show 
several possible realizations of this mechanism getting sufficiently stiff equations of state. The 
appearance of hyperons in dense neutron star interiors is accounted for. The obtained equations 
of state remain sufficiently stiff if the reduction of the <j) meson mass is incorporated. Thereby, 
the hyperon puzzle can be resolved. 


1. Introduction 

In this contribution we present the results of our studies of the equation of state (EoS) for the 
cold nuclear matter. We work on the extension of a relativistic mean-held (RMF) model [T] 
including a possibility for scaling of hadron masses and meson-baryon coupling constants. As 
in [T|, all hadron masses are scaled in the same way, whereas the coupling constants are scaled 
differently. 

The current challenge is to construct such an EoS that can describe a neutron star (NS) with 
the mass of (2.01 ± 0.04)Mq, where Mq is the mass of the Sun, which was recently observed, 
cf. [21 [3]. For that the EoS should be rather stiff. Simultaneously the EoS has to respect the 
constraint derived from the how of particles produced in heavy-ion collisions [ 1 ] . That oppositely 
requires the EoS to be rather soft. The additional complication in the description of heavy NSs 
arises if hyperons are included into consideration, since their presence softens the EoS, cf. |5]. 


2. The RMF model with scaling 

The nuclear energy-density functional of the RMF model including members of the baryon SU(3) 
octet {B = N, A, S, H) which interactions are mediated by the scalar meson (cr) and non-strange 
members of the vector-meson SU(3) nonet V = uj^p^cj), is 
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Here py,b = is the baryon Fermi momentum with ns standing for the density of 

baryon B. The effective baryon density fiv, being a source for the mean-field of meson V, 
is defined as ny = YIb ^VBQvnB, where qy stands for the baryon charge for H = w, the 
baryon third isospin projection for V = p, and the baryon strangeness for V = (f). The effective 
baryon mass m*^ differs from the vacuum mass niB because of the coupling to the scalar field, 
~ — XcrB^riN f 1 where / = gaN(x/mN is the dimensionless scalar field, and g^N is the 

(tN coupling constant. The functions gy{f) are the key ingredients, which differ our approach 
from other RMF models. These functions stem from the ratios of scalings of vector-meson masses 
to scalings of hadron-nucleon coupling constants and are tuned in our approach to render the 
best agreement of the resulting EoS with empirical constraints. The scalar field is found from 
the minimization of the functional ([T]) with respect to /: SE[f]/6f = 0. The coupling constants 
Ca,uip are to be determined from the saturation properties of the nuclear matter, therefore we 
have XaN = = XpN = 1- Since the (p meson does not couple to a nucleon, if the strict 

Okubo-Zweig-Iizuka selection rule is applied, i.e. = 0, the coupling constant C^j, is related 
to the coupling of the to meson as where = 783 MeV, = 1020 MeV. 

The coupling constants of vector mesons to hyperons are chosen traditionally according to the 
quark SU(6) symmetry 

XioA — Xi^j^ — — 3 , Xpj] — 2XpS — 2 , XfpA — XfpY, — , XpA — X(fi]\f — 0. (2) 

The coupling constants of hyperons with the scalar mean field are derived from the hyperon 
binding energies in isospin-symmetric matter (ISM) at the saturation density n = uq given 
by DP ^bfnd(”o) = C‘^xnj^XuiHno — mN+rn*j^{no) , and the empirical values = -28 MeV, 

^Wnd(^o) = 30 MeV and i’bind(”o) = —15 MeV. The repulsive S potential prevents the 
appearance of S hyperons in all models considered below. 

For the description of the NS one should add to m the energies of leptons. The NS 
composition is governed by the weak processes bringing the matter in the /3-equilibrium and 
relating chemical potential of all species to the electron and neutron chemical potentials. 
Together with the electroneutrality condition these relations allow to express the densities of all 
baryonic and leptonic species as functions of the total baryon density n = YIb The pressure 
of the NS matter follows then as P{n) = ndE/dn — E. With P = P{E) at hand, the NS masses 
and radii are calculated with the help of the Tolman-Oppenheimer-Voikoff (TOV) equation. In 
the crust region we match smoothly our EoS with the BPS EoS [6]. 

3. A novel mechanism how one can stiffen the EoS in an RMF model 

With appropriate choices of the scaling functions gy{f), and U{f) and, perhaps, also with a 
modification of our model can be reduced to various traditional RMF models. For 

example, putting 

gM) = gpif) = U if) = m%{bf/3 + cf^/P) (3) 

we recover the extension of the original Walecka model proposed by Boguta and Bodmer [3, 
the so-called non-linear Walecka (NEW) model, which allows to fit satisfactory the nuclear 
saturation properties, which we chose here for an illustration as follows: the saturation density 
no = 0.16fm“^, the binding energy per nucleon = —16MeV, the compressibility modulus 
K = 250 MeV, the effective nucleon mass m*j^{no)/mi\f = 1 —/o = 0.8, and the symmetry energy 
Esym = 30 MeV. For this choice of parameters and scalings ([3]) the limiting mass of a NS is 
1.92 Mq that is below the experimental limit. 

In [8] we proposed a simple method how one can stiffen the EoS in an RMF model. The 
strategy is to reduce the growth of the scalar field above a value /cut- As argued in [8], this can 
be achieved by adding to the original energy-density functional m a function of the scalar field. 



Figure 1. Results for the NLW and NLWcut models. Panel (a): the nucleon effective mass in 
the ISM as a function of the nucleon density. Panel (b): Pressure P as a function of the nucleon 
density in the ISM for various values Cg-; the insertion shows the binding energy per nucleon. 
The shadowed area is the constraint on the pressure from the particle flow in heavy-ion collisions 
(HIC) obtained in [1]. Panel (c): The NS mass as a function of the central density. The hatched 
band denotes the uncertainty in the value of the measured heaviest NS mass [3]. 


which is vanishingly small for / ~ /(uq) < /cut and is sharply increased for / ~ /cut) so that its 
second derivative becomes sufficiently large. Then, the self-consistent solution of the equation 
for / does not let the function /(n) to grow significantly above /cut within a broad interval of 
densities for n > n*, where n* is determined by /cut- We will illustrate this mechanism at hand 
of the following example: we replace 17 in ([3]) by U{f) + AU{f), where 

A17(/) = aln[l-Fexp(/3(/-/cut))] , a = ml, /3 = 120. (4) 

We will refer to such an extension as the “fj-cut scheme” and call the model as the NLWcut 
model. Applying the method to models with different values of /o it is convenient to parameterise 
/cut = /o + Co-(l — /o), and allow different values of to vary. The NLW model is recovered 
from the NLWcut model for Co- ^ 1 corresponding to /cut ~ 1- 

On panel (a) in Fig. [T]we show the effective nucleon mass m*j^ as a function of the nucleon 
density in ISM. We see that the function flattens off for densities n > 1.9no, 

3.0 no, and 4.0no for = 0.2, 0.3, and 0.4, respectively. The influence of the quenching of the 
nucleon-mass decrease on the EoS is illustrated on panel (b) in Fig. [H where we show the total 
pressure P{n) and the binding energy per nucleon as functions of n in the ISM. We observe that 
the replacement Q leads to a sizable increase of the pressure for n > n*. The shadow region 
shows the constraint on the pressure extracted from the particle flow in HIC [ 3 ] . We see that the 
EoSs in the NLWcut models satisfy this constraint for ^ 0.2. In difference from the pressure, 
the energy density and corresponding binding energy increase much weaker after the inclusion 
of the new term (j3]). The NS masses as functions of the central density are shown on panel (c) 
in Fig. [1] The stiffening of the EoS reflects in an increase of the NS mass at a given central 
density. The proposed cr-cut scheme allows us to shift the maximum NS mass from 1.92 Mq 
for the original NLW model to 1.96 Mq, 2.03 Mq, and 2.12 Mq for the NLWcut models with 
Co- = 0.4, 0.3, and 0.2, respectively. 

4. KVOR, KVORcut and MKVOR models 

In m the model KVOR was proposed with scaling functions rjp > 1 and < 1 for n > no, 
which were chosen to shift the threshold of direct Urea (DU) process, n ^ p + e + i?, to higher 
densities (via the r]p tuning) and to increase the maximum neutron star mass (via the tuning). 












Figure 2. Panel (a): The nucleon mass as a function of the density in KVOR, KVORcut and 
MKVOR models in BEM. Vertical bars indicate central densities in the stars with the maximum 
masses for the corresponding model. Panel (b): Pressure as a function of the nucleon density 
for ISM. Double-hatched area is the constraint from the particle flow in heavy-ion collisions [3]. 
Panel (c): The NS mass versus the central density. Bold dots show maximum available masses. 
The band shows uncertainty range of the measured mass of PSR J0348-I-0432 [3]. Horizontal 
dashes show the DU thresholds. 


The KVOR model exploits the saturation parameters m*j^{no)/m]\f = 0.805, K = 275 MeV, 
Jo = 32 MeV, and other parameters being the same as for the NLW model discussed above. The 
KVOR EoS was shown in [9] to satisfy appropriately the majority of experimental constraints 
known by that time. Hyperons were not included in the original KVOR model. However even 
without hyperons the KVOR EoS yields the maximum mass 2.01Mq that fits the constraint [3] 
only marginally. 

We discuss now several extensions of the KVOR model [T], which can better fulfill the modern 
constraints on the nuclear EoS. To stiffen the EoS we use the method described in the previous 
section but now will change not the potential U{f) but the function as 

^ + Y [1 + tanh(6.(/ - /,„,.))] . (5) 

We will call such extensions the KVORcut models. We consider three choices of parameters 
Qij, bco and /cut,a;: KVORcut04 with a^i = —0.5, = 55.76, /cut,a; = 0.454; KVORcut03 with 

= -0.5, b^ = 46.78, /cut,a; = 0.365 ; KVORcut02 with = -0.2, b^ = 74.55, = 0.249. 

In [1^ we introduced another RMF model with scaling — the MKVOR model — where the 
mechanism of the / stabilization is implemented not in the function but in the r]p(f) 

function, which rises for / < 0.5 similarly to the KVOR model but drops fast for / > 0.52. 
Therefore in the MKVOR model the nucleon mass decreases with the density increase in the 
ISM matter stronger than in the /3-equilibrium matter (BEM). All other scaling functions we 
chosen anew to satisfy better various constraints on the EoS form microscopic calculations, 
see m for details. 

The nucleon effective masses for KVORcut and MKVOR models are shown in Fig. [2l panel 
(a), as functions of the baryon density for BEM. In all models the effective hadron masses 
hrst decrease with the density and then saturate at some approximately constant values. The 
pressure as a function of the nucleon density in the ISM is shown on panel (b) in Fig. [2j We 
see that pressure for the KVORcut04 model is only slightly higher than that of the KVOR one. 
The curve for the KVORcut03 model goes on top of the upper boundary of the experimentally 
allowed region. The MKVOR-model curve goes through the region for n < 4no but escapes it 
at higher density. Oppositely, the curve for the KVORcut02 model lies above the region for 




















n < 3.5no but enters it for higher n. Panel (c) in Fig. [2] demonstrates the NS mass as a function 
of the central density for our EoSs. The band shows the uncertainty range of the masses for 
PSR J0348+0432 [3] which is fulfilled by all models. We see that for the KVORcut models the 
lesser the value is chosen, the larger the value of the maximum mass is and the smaller the 
central density corresponding to M^ax becomes. The MKVOR EoS is the stiffest among the 
EoSs considered here. Dashes on the right panel indicate the DU thresholds. The MKVOR 
model yields the highest DU threshold density. Thus, the KVORcut03 and MKVOR models 
are the most promising models for the simultaneous fulfillment of the particle-flow and NS-mass 
constraints. 

5. Inclusion of hyperons 

We discuss now the results of the hyperonization phase transition on the EoS within our RMF 
models. We will consider three schemes for the inclusion of hyperons: 

(i) The hyperons are included with the vector-meson coupling parameters XujH and XpH according 
to m and the scalar coupling parameters x^h constrained by the hyperon binding energies as 
outlined in section [2j We switch off the hyperon coupling to cj) mesons, putting x^pn = 0. The 
models with hyperons obtained within such a scheme are labeled as the KVORH, KVORHcut03 
and MKVORH models. 

(ii) Next we incorporate the (/>-meson mean field and assume the very same scaling of its mass 

as for all other hadrons, i.e., rn*^{f) = — /), but with the unsealed hyperon-^ coupling 

constants, which gives the scaling function //</,(/) = (1 — /)^. The resulting models are denoted 
as KVORH(/), KVORHcut03(/) and MKVORH(/). 

(hi) In addition to the (^-meson mass scaling we take into account a possible change of the 
XaH parameters with an increase of density. We assume these parameters being unchanged 
for densities n ^ no and decreasing for higher densities, so that they reach and stay zero for 
the baryon density n > UcH, where UcH are critical densities for hyperonization. With such a 
parametrization we will exploit vacuum masses of the hyperon H for n > UcH- Such extensions 
of the models are call KVORHcut03(^cT and MKV0RH(()(T models. 

As known, the inclusion of hyperons results usually in a substantial decrease of the maximum 
NS mass. The difference between NS masses with and without hyperons proves to be so large for 
reasonable hyperon fractions in the standard RMF approach that in order to solve the puzzle one 
needs to start with a very stiff EoS without hyperons that would be in odds with microscopical 
calculations and with the constraints from heavy-ion collisions. 

The NS masses as functions of the central baryon density for various models with hyperons 
are shown in Fig. [3j Comparing with Fig. [2] we see that the appearance of hyperons reduces 
the maximum NS mass for the KVORH and KVORHcut03 models by 0.35 Mq and 0.33 Mq, 
respectively, bringing them below the empirically known pulsar mass. We find also that hyperons 
cannot be incorporated in the MKVOR model, since the solution for f{n) does not exist for all 
densities. Including the (/>-meson with an /-dependent mass and scaling of aH coupling constants 
(schemes (ii) and (iii) above) we obtain maximum NS mass of 1.97 Mq for KVORHcut03/ that 
marginally agrees with the observational constraint and 2.07 Mq for KVORHcut03())(T that fully 
fulhls the maximum NS mass constraint. For the extensions of the MKVOR models we find the 
maximum mass of 2.22 Mq for MKVORH/) and 2.29 Mq for MKVORH/cr; both values are well 
above the observational maximum NS mass constraint. 

Mass-radius relations for KVORHcut03/), KVORHcut03(/)(T, MKVORH/ and MKVORH/cr 
EoSs are shown on panel (c) in Fig. [3l We see that the MKVORH/ and MKVORH/a EoSs 
satisfy all presented constraints, whereas the KVORHcut03/ and KVORHcut03/(T EoSs render 
star radii slightly larger than the corridor of allowed values (shown by black thin lines) extracted 
in m from the analysis of empirical data. 



Figure 3. Panels (a) and (b): NS mass versus the central baryon density for various models 
with hyperons. Panel (c): Mass-radius relation for KVORcutOS and MKVOR models including 
hyperons and constraints from thermal radiation of the isolated NS RX J1856 m and from 
quasi-periodic oscillations in the LMXBs 4U 0614-1-09 [12]. The black thin lines show the radius 
constraint from the analysis m- The horizontal band shows the uncertainty range of masses 
for PSR J0348-h0432 |3]. 


6. Conclusion 

We presented a novel mechanism of stiffening the EoS within an RMF model with hadron masses 
and coupling constants depending on the scalar field. The key idea is to quench the growth of 
the scalar field and correspondingly a decrease of the nucleon mass with a density increase. 
Three different approaches are outlined where the quenching is achieved by the modification of 
the a mean-field self-interaction or by changing in the vector (w) and isovector (p) sectors of the 
energy density functional. Within these schemes we discuss the cr-cut versions of the traditional 
non-linear Walecka model, the ca-cut versions of the KVOR model proposed in [T] (KVORcut 
models) and a new model MKVOR combining various schemes. The hyperons are included and 
it is demonstrated that KVORcut and MKVOR models can support sufficiently heavy neutron 
stars if a decrease of the (/>-meson mass or the hyperon-cr coupling constants are incorporated. 
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